Quantum phase transitions in the systems of parallel quantum dots 
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We study the low-temperature transport properties of the systems of parallel quantum dots described by the 
iV-impurity Anderson model. We calculate the quasiparticle scattering phase shifts, spectral functions and 
correlations as a function of the gate voltage for N up to 5. For any N, the conductance at the particle-hole 
symmetric point is unitary. For N > 2, a transition from ferromagnetic to antiferromagnetic impurity spin 
correlations occurs at some gate voltage. For N > 3, there is an additional transition due to an abrupt change 
in average impurity occupancy. For odd TV, the conductance is discontinuous through both quantum phase 
transitions, while for even TV only the magnetic transition affects the conductance. Similar effects should be 
experimentally observable in the systems of quantum dots with ferromagnetic conduction-band-mediated inter- 
dot exchange interactions. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.63.Kv 
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Parallel double quantum dots in semiconductor heterostruc- 
tures exhibit many interesting quantum effects at low temper- 
atures, such as the formation of molecular states, Aharonov- 
Bohm oscillations, phase lapses and the Kondo effect |H 0, 
01 . They are also predicted to exhibit quan- 
tum phase transitions of different kinds I0, |9[]. Inter-dot 
exchange interactions, both ferromagnetic (FM) and anti- 
ferromagnetic (ATM), play a central role in such systems 
lES [Ull [H SI Gj, M O. Furthermore, the phase co- 
herent electron transport leads to various possible interference 
effects and Fano anti -resonances ill 8l Il9l l20l 12111 . Recently, 
few-electron triple quantum dot structures have been fabri- 
cated 11221 12311 and even more complex multi-dot nanostruc- 
tures can in principle also be assembled. It is thus appropri- 
ate to endeavor the studies of new phenomena which occur in 
multiple-dot systems. 

The simplest A^-impurity Anderson model for several dots 
embedded in parallel between two conduction leads in a left- 
right symmetric way (see the insets in Fig. d) is defined by 
H = Hb + Yli=i ^i, where describes a band with a con- 
stant density of states p = 1/2D (2D is the bandwidth) and 



Hi = 5rii 
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H.c). (1) 



Parameter 6 = e + U/2 is related to the gate voltage, U is 
the electron-electron repulsion and we assume that all impu- 
rities hybridize with the same left-right symmetric combina- 
tions of states from both leads with a constant hybridization 
function T = irpV 2 ; this corresponds to taking the limit of 
small inter-dot separation [8]. The inter-dot tunneling cou- 
pling and capacitive coupling (inter-dot charge repulsion) are 
assumed small and all dots are equivalent: the system thus 
has Sn symmetric group symmetry of all possible permuta- 
tions of dot labels i. At the particle-hole (p-h) symmetric 
point, 6 — 0, and for U/ttT 3> 1, the conduction-band- 
mediated inter-dot exchange interaction induces FM align- 
ment of the impurity spins and the system undergoes the 
spin-A^/2 Kondo effect ending up in an underscreened strong- 



coupling (SC) Fermi liquid fixed poin t with residual spin 
N/2 - 1/2 |i EE [H El H SlIS 0. For very large 
5/U, the impurities are unoccupied and the system is in the 
frozen-impurity (FI) fixed-point with no residual spin. In the 
single-impurity (N = 1) case, the SC and FI fixed points lie 
on the same line of fixed points and they differ only in the 
strength of the potential scattering [2f|]. For N > 2, however, 
the SC and FI lines of fixed points are qualitatively different 
(each corresponding to a different residual spin) and must be 
separated by at least one quantum phase transition (QPT) B29I1 . 

A reliable technique to study the low-temperature proper- 
ties of coupled quantum dot systems is the numerical renor- 
malization group (NRG) [30, 31]. In this Letter, we show that 
for any N > 2, the FM alignment collapses at some critical 
value S c i and that for 5 > 5 C \ the inter-impurity spin-spin 
correlations are AFM. For > 3, there is precisely one addi- 
tional QPT at slightly higher 5 C 2 related to an abrupt change in 
the average impurity occupancy. The essential novelty of our 
results is that certain of these phase transitions can be easily 
detected in zero-bias conductance measurements. 

Conductance. - The on-site energy e can be regulated using 
the gate voltages to tune the charge state (occupancy) on the 
dots. Gate-voltage dependent conductance is shown in Fig.QJ 
for N — 1, . . . , 5 for a range of magnetic field strengths. The 
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conductance is calculated as G = Gq/2^2^^ 
where Go = 2e 2 /h is the conductance quantum 1 32] and the 
quasiparticle scattering phase shifts 5° p are extracted from 
the NRG excitation spectra. 

At S = 0, the systems are fully conducting at zero field 
and there is a wide plateau of high conductance associated 
with the spin-AT/2 Kondo effect fl. While the A^ = 1 sys- 
tem smoothly crosses over from the Kondo regime to the non- 
conducting FI regime, in the multi-impurity case we observe 
sharp discontinuities: one discontinuity for A^ even and two 
discontinuities for A^ odd. The conductance culminates in a 
unitary peak slightly below e = (i.e. below 5/U = 1/2) for 
all N > 2. The origin of this peak is simply potential scat- 
tering. The magnetic field B has a strong effect on the Kondo 
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Figure 1: (color online) a) Zero-temperature conductance through 
systems of N parallel quantum dots as a function of the gate volt- 
age for a range of magnetic fields. The in-plane field only leads to 
Zeeman splitting; no magnetic flux pierces the rings formed by pairs 
of dots. A = 4; NRG iterations were performed until the zero tem- 
perature limit was reached. T K = 3.3 1(T 6 D (for N = 2). The 
magnetic field is measured in units of g^B- Only S > is shown due 
to the symmetry of the problem, b) Zero-temperature phase diagram 
delimiting the different regimes as a function of the gate voltage. 
Filled circles (•) correspond to phase transitions visible in conduc- 
tance, while the empty circle (o) denotes the phase transition with 
no associated conductance discontinuity, c) Conductance and spin 
conductance in small magnetic field in the transition region. 



plateau: the conductance is significantly reduced as soon as 
B is of the order of the Kondo temperature Tjc- The potential 
scattering peak is only affected by extremely high fields of the 
order of U. 

Quantum phase transitions. - Conductance discontinuities 



find their counterparts in the jumps of the total impurity oc- 
cupancy and spin-spin correlation (Si ■ Sj) (Fig. [2^); a new 
feature, however, is the existence of two points of discontinu- 
ity for N = 4 while the conductance only exhibits one. In the 
Kondo regime for 6 < <Li, the systems are nearly half-filled 
and spins are aligned [8]. As we cross 5 C \, the occupancy 
abruptly decreases and the spin correlations turn from FM to 
AFM. For A > 3, a second discontinuity occurs at somewhat 
higher 5 C 2\ its characteristic property is the occupancy jump 
by almost exactly N — 2, from N — 1 to 1. According to the 
Friedel sum rule, a change in the occupancy by n is mirrored 
in a change of the scattering phase shift by A<5 q . p . = nir/2. 
This explains the conductance jump from G = to G = Go 
in the case of odd N > 3 and the absence of the second con- 
ductance discontinuity for even N > 4. It is remarkable that 
the second QPT occurs precisely at the point where the con- 
ductance is extremal. 



a ) U/D=l, I7U=0.04 
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Figure 2: (color online) a) Total occupancy (charge) and spin-spin 
correlation between pairs of spins as a function of the gate voltage, 
b) Schematic representations of the spin configurations in the inter- 
mediate AFM ordered phases, their total spins and degeneracies for 

N > 3. 



The discontinuities originate from an interplay of the 
RKKY interactions [^L the presence of the bound states in 



the continuum 1JT 




and the occupancy switching 



sins. We note that only the symmetric state described 
by the operator d| ym = 1 /VN J^i hybridizes directly with 
the conduction band, while all asymmetric states are decou- 
pled. A close-up on the points of discontinuity for N = 2, 3 
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is shown in Fig. [3] As we move away from S — 0, only the 
symmetric state is depopulated until we reach At this 
point, the asymmetric levels are depopulated and the symmet- 
ric level re-populated (up to exactly 1 for A = 2); such oc- 
cupancy switching stems from the competition between the 
T-dependent coupling of dots to the conduction band and the 
charging energy U B35I1 . Between S c \ and 5 C 2, the occupancy 
of all states decreases until at 5 C 2 another charge oscillation 
occurs in which the occupancy of the asymmetric states plum- 
mets. Both transitions can be classified as first-order boundary 
quantum phase transitions (level-crossing) 03311 : the schematic 
phase diagram in shown in Fig.QJ. 



in the continuum) found in the systems of noninteracting par- 
allel dots mim. 
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Figure 3: (color online) Expectation values in the symmetric and 
asymmetric basis for N = 2 and 3 systems. The asymmetric combi- 
nation of states considered is dJ. S y m = 1/V2(d\ — d\). 

For S cl < 6 (A = 2) and for S cl < S < 5 c2 (A > 3), 
each spin interacts with every other spin antiferromagneti- 
cally with an equal strength. The ground state of such ef- 
fective Heisenberg Hamiltonian consists of one singlet for 
A = 2, two degenerate doublets for N = 3, two degener- 
ate singlets for N = 4, etc. [36], see Fig. |2j). Results for 
magnetic susceptibility and entropy demonstrate that in odd- 
A systems the spin degree of freedom is screened by the spin- 
1/2 Kondo effect and that the residual entropy is the logarithm 
of the additional degeneracy of the ground state spin multi- 
plets. This implies that in this phase (for A > 3) there are 
several equally probable ways for the electron spin ordering. 
Consequently, this "magnetic-frustration" phase is sensitive 
to the breaking of the Sn symmetry between the impurities 
which lifts the degeneracy. For odd A and for weak magnetic 
field, the transmitted current in this regime becomes fully spin 
polarized at some gate voltage (it should be noted that in the 
presence of magnetic field the phase transition at 5 C 2 is re- 
placed by a cross-over 12711 '). see the plot of G sp i n /G where 



Gspin = Go/2(sin 2 — sin 2 S^ p ) is spin current and G is 
total (charge) current, Fig, lfc. Thus the A — 3 device might 
function as a spin filter ||27 ^37 , 381. 

For decreasing interaction strength U, the discontinuous 
features in the conduction plots narrow down and disappear 
for U — 0; at this point the conductance curve consists of a 
single Lorentzian peak of width AT centered at 5 = 0. This 
demonstrates that the conductance discontinuities are adiabat- 
ically connected with the ghost Fano resonances (bound states 



Robustness. - Finite inter-dot charge repulsion U\2 and 
inter-dot hopping t have little effect as long as U12 <C U, 
t -C U - U12 and At 2 /(U - f7 12 ) < J e g, where J cff is the 
effective conduction-band mediated inter-impurity exchange 
interaction [8]. The first two constraints are easily met in ex- 
periments fl39ll and together they imply the third unless J c ff is 
very small. 

The Sn symmetry is broken if the strengths of the hy- 
bridization Ti of each impurity to the conduction band are 
made unequal, or if different gate voltages are applied to 
the dots. This leads to smearing of the discontinuities in the 
occupancy and correlation functions; nevertheless, the discon- 
tinuities in the conductance curves persist (i.e. level-crossings 
still occur). Curiously, the first discontinuity in the conduc- 
tance no longer coincides with the sign change of (Si ■ Sj). 

For a generic problem of parallel dots, in particular if left- 
right (L-R) and Sn symmetries are weakly broken, the anti- 
symmetric combinations of the conduction band electrons be- 
come relevant and we need to consider the full A-impurity 
two-channel model. In this case, the spin- A/2 Kondo ef- 
fect is followed by another stage of the Kondo screening to 

(2) 

S = A/ 2 — 1 at significantly lower Kondo temperature Ti- . 
This leads to a phase shift of 7r/2 in the odd scattering channel 
and the zero-temperature conductance in the "Kondo regime" 
becomes very small Q (see Refs. H ESS for the A = 2 
case). Nevertheless, if L-R breaking is weak, the conduc- 
tance curves shown in Fig. Q] are a good approximation for the 
finite-temperature conductance in the experimentally relevant 
T K > T > Tg } range. 

Spectral functions. - In Fig.Hk, we plot spectral functions 
A lj (u) = -l/(27r)Im(G^ + G^) for N = 1 and 2 at the 
p-h symmetric point |40]. The diagonal spectral functions 
Aj = An represent the on-site density of states, while the 
out-of-diagonal spectral functions Aij with i 7^ j are related 
to the processes where one electron is injected at one site and 
later extracted at a different site. The peak in Ai at U/2 is the 
familiar charge excitation peak that is observed for all A; for 
A > 2 a small negative peak appears in the out-of-diagonal 
spectral densities Aij . Additional features for J e g < u> < U 
are related to the magnetic alignment the diagonal 

spectral function exhibits a broad hump which peaks at ~ J c g, 
while the out-of-diagonal spectral function exhibit a slight de- 
pression. For u> < J c s when the spins align, all A^ curves 
merge into a Kondo resonance. 

In Fig. UK we plot the symmetrized and normalized spectral 
density function g(ui) = Aij, which determines the 

conductance through the dots as G = Go J (—d f/doj)g duj, 
where / is the Fermi-Dirac distribution function fSBl. In 
a simple approximation, g(u>) at u> = T indicates the conduc- 
tance through the system at temperature T. Only for A = 1 
is the approach to the unitary conductance at T = rapid 
(quadratic), as expected for regular Fermi liquid systems. For 
A > 2, the Kondo resonance is cusp-like and the approach 
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Figure 4: (color online) a) Normalized spectral functions nVAij(u>) 
and b) function g(ui) = 7iT }[V . j4ij(oj). 



to the unitary conductance is very slow (logarithmic) [ 24, 43]. 
This is characteristic for underscreened Kondo systems which 
behave as singular Fermi liquids ll43il44ll . 

Conclusion. - For N > 3, the A^-impurity Anderson model 
undergoes two phase transitions. The first transition separates 
the spin alignment and the associated spin- N/2 Kondo screen- 
ing from the spin anti-alignment with magnetic frustration and 
(for odd N) Kondo screening of the spin- 1/2 moment. The 
second transition reflects the instability of the phases with the 
occupancy in the interval 1 < (ntot) < N — 1. Furthermore, 
for odd N the system abruptly switches from being fully con- 
ducting to zero conductance; this would facilitate the experi- 
mental observation of similar effects in quantum dot systems 
and might even have applications as a switch or a transistor 
with very high on-off ratio. In addition, the odd A" system ap- 
pears as a possible realization of gate voltage switchable spin 
filter device. 
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